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Electronic Raman scattering in high-T c superconductors is studied within the t — J model. It is 
shown that the Ai g and Bi g spectra are dominated by amplitude fluctuations of the superconducting 
and the d-wave CDW order parameters, respectively. The £?2 9 spectrum contains no collective effects 
and its broad peak reflects vaguely the doping dependence of T c , similarly to the pronounced peak 
in the A\ g spectrum. The agreement of our theory with the experiment supports the picture of two 
different, competing order parameters in the underdoped regime of high-T c superconductors. 
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Electronic Raman scattering in superconductors 
probes charge excitations across the superconducting gap 
and thus provides information on the magnitude and 
the anisotropy of the gap [1—4]. Commonly used weak- 
coupling theories for the interpretation of the data are, 
however, not able to account for the experimental spectra 
in the cuprates, especially, as a function of doping. Such 
theories yield values for the d-wave gap which have no 
simple relation to T c and which are substantially smaller 
than those measured in ARPES experiments [5]. In our 
view, these problems arise because electronic correlation 
effects are not adequately taken into account. Below we 
consider a strong-coupling model which allows to con- 
sider interaction effects between excited quasi-particlcs 
as well as the appearance of a d-CDW at lower dopings 
in a systematic way. 

Our calculations of Raman response functions are 
based on the widely accepted t-J model where the two 
spin components have been generalized to N compo- 
nents and the leading diagrams at large N are taken 
into account. As discussed in detail in Ref. [6] the 
phase diagram in this limit is largely determined by 
the onset of a flux phase with d-wave symmetry, of- 
ten also called d-CDW [7], at a doping 5 = So- The 
corresponding order parameter has the form <fr(k) = 
-i/2A c J] q(T J(k - q)(c qCT c q+QCT ). J is the Heisenberg 
coupling, , c are creation and annihilation operators 
for electrons under the constraint that double occupan- 
cies of lattice sites are excluded, N c is the number of 
primitive cells, (...) denotes an expectation value, and 
Q is the wave vector of the d-CDW. Furthermore, there 
exists for all dopings an instability towards d-wave su- 
perconductivity [8]. Keeping only the instantaneous 
term in the effective interaction, the order parameter 
is A(k) = l/2iV c £ (J(k - q) - V c (k - q))(£ qT 5_ ql ). 
As shown in Ref. [6] it is in general necessary to in- 
clude the Coulomb potential Vc in order to stabilize the 
d-CDW with respect to phase separation. The result- 
ing superconducting transition temperature T c decreases 



with increasing doping for 6 > Sq. For S < Sq the two 
order parameters compete with each other leading to a 
strongly decreasing T c with decreasing doping. As a re- 
sult optimal doping is essentially determined by Sq. In 
the presence of the two order parameters the operators 
(c^ *, C-kjj, Cjj^q *, c_k-Q,j,) are coupled leading to the 
following Green's function matrix [6] 

/z-e(k) -A(k) -i$(k) 0_ \ 

1( -A(k) z + e(k) i$(k) 

^> K J- i$(k) 0_ z-e(_k) -A(k) 

\ -t$(k) -A(k) z + e(k) / 

(1) 

e(k) is the one-particle energy, e(k) = — (St + 
aJ)(cos(k x ) + cos(ky)) — 2t' Scos(k x )cos(k y ) — fx, with 
q = 1/N C C0S (Qx)I{ e ( ( l))- / is the Fermi function, 5 
the doping away fom half-filling, fi a renormalized chem- 
ical potential, t and t' are nearest and second-nearest 
neighbor hopping amplitudes, z a complex frequency, and 
k = k Q. 

Expressing the expectation values in the order param- 
eters by G and using Eq.(l) one obtains two coupled 
equations for the order parameters. Detailed consider- 
ations [6] show that for dopings in the neighborhood of 
the optimal doping the most stable order parameters 
have d-wave symmetry, <I > (k) = $7(k),A(k) = A7(k), 
with 7(k) = (cos(k x ) — cos(k y ))/2. Fig. 1 shows the 
doping dependence of $ and A at zero temperature, cal- 
culated for t'/t = —0.35 and J/t = 0.3. The energy unit 
is t. A repulsive nearest-neighbor Coulomb interaction 
was also included with Vc/t — 0.06. In the overdoped 
region 5 > Sq $ is zero and A increases monotonically 
with decreasing 5. After the onset of the flux phase at 
Sq $ is first suppressed by superconductivity and then, 
with decreasing doping, increases steeply and suppresses 
now, in addition to V, the superconducting order param- 
eter. In our calculation we assumed the value (it, tt) for 
the wave vector Q. Fig. 1 shows that superconductivity 
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coexists in the underdoped region with the d-CDW and 
that the competition between the two order parameters 
leads to the rapid decay of T c towards small dopings. 
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FIG. 1. Order parameters $ and A as a function of doping 
in units of t at T = 0. 

The scattering intensity of electronic Raman scattering 
is proportional to — (1 + n(w))/mx a (w + it]), where n 
denotes the Bose factor, u) the frequency and \ a a re- 
sponse function of the non-local density operator p a (k) = 
1 /^cI] k . (T 7a( k )c^ (T 5k<T- a = 1,3,4 stands for the rep- 
resentations Ai g , Big, B 2g of the point group D±h of the 
square lattice, respectively. In the large N limit of the 
t-J model we have 71 (k) = (St + Ja)(cosk x + cosk y )/2 + 
2St'cosk x ■ cosky, 73(k) = (St + Ja)(cosk x — cosk y )/2, 
74 (k) = —2St' sink x ■ sink y . The response function \a is 
determined at large N by the sum over ladder diagrams 
shown in Fig. 2. 




FIG. 2. Ladder diagrams for the Raman susceptibility; the 
filled triangle is the Raman vertex j a , the solid and dashed 
lines represent the Green's function G and a residual interac- 
tion, respectively. 

In the case of the B lg and B 2g spectra the residual 
interaction is the Heisenberg interaction. Decompos- 
ing it into irreducible basis functions of D ih it contains 
one Big kernel but no contribution with B 2g symme- 
try This means that \i is given by the free susceptibil- 
ity which reduces to the standard expression if the 
d-CDW order parameter vanishes. We find the follow- 
ing analytic expressions x 3 (» = xf\ z )/{ 1 + J'xf\z)), 
J' = J/(St + Ja) 2 , X4 (z) = X i 0) (z), with 

x« } W = w £ 7 «W ( n iMi( k > z) n 12;21 (k, z) 



+ (-l) a [n 13 , 3 i(k, z) - n 14 ,4i(k, «)]) , (2) 
U ijt ki(Kz) = T~y]Gij(iu n + z,k)Gki(ioJ n ,k), (3) 

n 

for a = 3,4. 

Fig. 3 shows Bi g and B 2g spectra for the parameters 
t' It = -0.35, J jt = 0.3, and T = 0. The upper dia- 
gram, corresponding to S = 0.178, A = 0.040 and $ = 0, 
is typical for the overdoped region. In the B 2g spectra 
74 heavily weights transitions near the diagonals in k- 
space where the gap vanishes. As a result the spectrum is 
broad and peaks substantially below 2 A. Neglecting ver- 
tex corrections x^ is determined by free particle-hole ex- 
citations mainly across the maximum of the d-wave gap. 
This leads to a well-pronounced peak at 2A as shown 
by the dashed curve in the upper diagram of Fig. 3. In- 
cluding also vertex corrections X 3 develops a bound state 
well within the d-wave gap which is somewhat broadened 
by the finite density of states in the d-wave gap. Practi- 
cally all spectral weight is shifted from the region around 
and above 2A into the bound state which, in anlogy to 
semiconductors, can be viewed as an exciton state [9]. 

The middle and lower diagrams in Fig. 3 decribe a 
slightly and a strongly underdoped case, respectively. In 
the extreme case of $ > > A the density of quasi-particles 
shows a small gap due to superconductivity embedded 
into a larger gap structure due to the flux phase. The 
B 2g spectra (dotted lines) are always very broad with 
a maximum which is located roughly in the middle of 
the BCS-gap which decreases with decreasing doping. 
In addition it contains a weak and structureless back- 
ground which extends over the whole gap region. The 
dashed lines in Fig. 3 describe the B^J spectrum due 
to non-interacting particle-hole excitations. It is rather 
insensitive to the BCS-part of the gap and consists of 
one strong peak around the maximal total gap. Taking 
also vertex corrections into account practically the total 
spectral weig ht of the B[ J curve is shifted into one peak 
which monotonically increases with decreasing doping. 
It describes amplitude fluctuations of the d-CDW order 
parameter. 

The essential features of Fig. 3 agree with experiments 
in the cuprates, in particular, the increase of the fre- 
quency of the Big peak with decreasing doping and the 
nonmonotonic behavior of the B 2g peak as a function of 
doping similar to that of T c [3]. Quantitative fits to ex- 
perimental curves [10] need, however, the inclusion of ad- 
ditional interactions, e.,g., with impurities, and are there- 
fore not attempted here. To what extent the predicted 
large excitonic effect in the B ig spectrum is compati- 
ble with experiment is presently unclear: The larger gap 
values deduced from ARPES data in comparison to the 
Raman Bi g peak in the 
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FIG. 3. Electronic Raman spectra of B-symmetry for three 
different dopings S, calculated for t' /t = —0.35, J/t = 0.3, and 
T = 0. The B29 spectrum has been multiplied by 10 in the 
two upper and by 100 in the lower diagrams. 

undcrdopcd regime and the small difference between B\ g 
and Big peaks on the overdoped side are in favor of exci- 
tonic effects in the B\ g channel. ARPES gives in slightly 
overdoped £>z2212 2A ~ 70 — 75meV [5] whereas the 
Big peak lies near 59 — 63meV [1,4], for the same dop- 
ing. ARPES yields in untwinnned Y BaiCuO'j gaps of 
58 and 88 meV at the X and y-points, respectively [11], 
wheras the corresponding peaks in Raman scattering arc 
at 50 and 60 meV [12]. Excitonic effects in the B\ g chan- 
nel thus seem not to be in conflict with available data. 
The overall decrease in intensity as a function of doping 
of the theoretical B-spectra is mainly caused by the 5 2 
factor in the Raman vertices. It agrees with experiment, 
however, the experimental B\ g spectrum decreases faster 
than the B^g spectrum in contrast with the theoretical 
curves. 

The dotted lines in Fig. 4, called Ai giS , show screened 
Ai g spectra. Their absolute intensities are about 2 or- 
ders in magnitude smaller than those for B± g scatter- 
ing. Moreover, they show in general two gap features, a 
weaker one related to the BCS and a stronger one reflect- 
ing the total gap. All these features clearly indicate that 
the Ai gtS curves are unable to explain the experimental 
Ai g spectra. 

According to Ref. [13] density fluctuations may couple 
to the superconducing order parameter via the modula- 
tion of the density of states at the Fermi energy. This 



indirect coupling also applies to the density fluctuations 
induced by the the ^4i g Raman vertex 71. Using the 
BCS-assumption with a cutoff uj the effective Raman 
coupling is J2k 7i(k)CkT^-ki + h- c -; with 71 = <77(k) and 
g = Aln(w /A)6N(0)/N(0), where 6N(0) is the change 
in the projected d-wave density at the Fermi surface due 
to the original Raman vertex 71. The effective vertex 
71 is unscreened. Moreover, g is large in the optimal 
and overdoped region because of the proximity of the 
van Hove singularity which enhances SN(0) and decays 
rapidly in the underdoped region because of the factor 
A. The resulting A\ g spectrum is again given by the 
diagram in Fig. 2 with external Raman vertices 71 and 
the Heisenberg interaction for the dashed line. Denot- 
ing the resulting susceptibilities with a tilde we obtain 

*1 0) W = W E^( k )( n !242(k, Z) + n n ,22(k, Z ) 

JVC k 

-IIi4,32(k, z) - n 13 ,42(k, z) + (z -» -z)). (4) 

The dashed and solid lines in Fig. 4, denoted by A^ g 

and Ai g , show the negative imaginary part of Xi°^ an d, 
Xi, respectively, for the same three dopings as in Fig. 3, 
calculated with the cutoff lo = J. A^j exhibits a step- 
like feature at 2A in the overdoped regime which trans- 
forms into a well-pronounced peak towards lower dop- 
ings. The peak coincides with the maximum of the total 
gap and there is no indication of any contribution from 
the smaller BCS gap, similar as in the B^J spectrum. 
Including also vertex corrections spectral weight around 
2A accumulates into a pronounced collective peak at 2A 
in the overdoped region describing amplitude fluctuations 
of the superconducting order parameter [13]. Going from 
the overdoped to the underdoped regime this peak first 
increases, then passes through a maximum around op- 
timal doping, and then decreases in frequency and in- 
tensity, becoming at the same time rather broad. The 
decrease in its frequency can easily be understood: The 
denominator of Xi (z) becomes zero at zero frequency in 
the limit A — > if A satisfies the gap equation. This 
implies that the frequency of the collective peak in Xi( z ) 
has to go to zero at the onset of superconductivity in 
the underdoped region in the absence of damping. The 
density of states is, however, nonzero in the d-CDW gap 
even at low frequency. As a result, the collective peak 
varies as 2A at larger energies but becomes rather broad 
at low energies. This agrees with experiments in Bi2212 
where the A\ g peak passes through a maximum and then 
decreases on the underdoped side [1]. Taking also data 
on phonon renormalizations into account it even has been 
conjectured [2] in the case of Y BaiCu-zOi that the A\ g 
peak corresponds to the superconducting gap 2A. 
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FIG. 4. A^J and Ai g : Unperturbed and full A\ g spectra 
calculated with the indirect coupling; Ai g y. Screened, full 
Ai g spectrum multiplied by 50 (upper two diagrams) or 10 
(lower digram). 

Fig. 4 shows that A\ g is much larger in magnitude 
than Aig :S for all dopings, in spite of the fact that the 
coupling to the superconducting order parameter occurs 
in an indirect way. One reason is that the 71 vertex is un- 
screened in contrast to 71 and that screening reduces the 
intensity by one order of magnitude or more. A second 
reason is that the proximity to the van Hove singularity 
leads to rather large values for 6N(0)/N(0) between 2 
and 5 in the considered doping region. Comparing the 
theoretical intensities for B\ g and A\ g scattering with 
the experimental ones one finds that A lg has the right 
magnitude whereas A lgtS would be much too small. We 
would like to point out, however, that the coupling con- 
stant g also depends on the cutoff cuq and thus is subject 
to considerable uncertainty. 

In conclusion, we have shown that the observed differ- 
ent behavior of the three symmetry components of the 
electronic Raman spectrum in high-T c superconductors 



as a function of doping can be explained within a t-J 
model in the larg N limit. Basic ingredients of this ap- 
proach are the strong competition of the superconduct- 
ing and the d-CDW order parameters in the underdoped 
regime and the importance of collective effects. The peak 
in the B lg spectrum in the superconducting state is ex- 
plained by amplitude fluctuations of the d-CDW order 
parameter which, in the optimal and overdoped region, 
can also be viewed as excitonic states. We also found that 
the indirect coupling of light to the superconducting or- 
der parameter is important leading to the conclusion that 
the Ai g peak is caused by amplitude fluctuations of the 
superconducting order parameter. 
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